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Abstract. We first summarize the classical arguments that the vast majority 
of glass-forming liquids require more than one "order" parameter for their 
description. Critiques against this conventional wisdom are then presented, 
and it is argued that the matter deserves to be reconsidered in light of recent 
experimental developments. Out of the eight basic thermoviscoelastic frequency- 
dependent response functions, there are generally three independent functions. 
For stochastic dynamics we show that there are only two independent response 
functions; for this case it is shown how analytic continuation may be utilized to 
express the third response functions in terms of two others. Operational criteria 
are presented for the linear thermoviscoelasticity being described by a single 
"order" parameter, in which case there is just one independent thermoviscoelastic 
response function. It is shown that a single "order" parameter description 
applies to a good approximation whenever thermal equilibrium fluctuations of 
fundamental variables like energy and pressure are strongly correlated. Results 
from computer simulations showing that this is the case for a number of simple 
glass-forming liquids, as well as a few exceptions, are briefly presented. Finally, we 
discuss a new conjecture according to which experiments at varying temperature 
and pressure follow the density-scaling expression for the relaxation time, t = 
F(p x /T) (p and T are density and temperature), if and only if the liquid is 
"strongly correlating," i.e., to a good approximation is described by a single 
"order" parameter. 
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1. Introduction 



The question whether one "order" parameter is sufficient for describing glass structure 
attracted considerable interest among the limited number of glass scientists in the 
period 1950-1980. The question was thoroughly discussed in particular in the 1970's 
[U0G2IU 5, 6 leading to clarifications of a number of theoretical questions. Since then, 
based on experimental evidence the consensus has been that one "order" parameter 
is rarely enough. 

The term "order parameter" was commonly used in the glass community before 
the term in the 1960's became commonly known in the physics community where it 
took on a somewhat different meaning. In connection with critical phenomena and 
the theory of second order phase transitions, renormalization, etc, "order parameters" 
reflect the relevant Lie group symmetry and determine the relevant part of the free 
energy within a standard Ginzburg-Landau expansion of the free energy. In order not 
to confuse the issue it is probably a good idea to change the wording, so below we 
refer to "order" parameters or occasionally just parameters. 

The present paper summarizes and extends recent works making the case 
that the question of how many "order" parameters are sufficient deserves to be 
reconsidered. In Sec. 2 we briefly summarize the classical viewpoint, in Sec. 3 
critiques against it are presented, in Sec. 4 the more restricted and well-defined 
case of linear thermoviscoelasticity is presented, in Sec. 5 thermoviscoelasticity in 
complete generality is discussed, in Sec. 6 we show that in any stochastic description 
of the dynamics there are only two independent response functions, Sec. 7 treats the 
single-parameter case where there is just one independent thermoviscoelastic response 
function, in Sec. 8 the new concept of a "dynamic" Prigogine-Defay ratio, which tests 
one-parameter-ness by reference to single-frequency thermoviscoelastic measurements, 
is presented. Section 9 presents a few computer simulations showing that several 
systems indeed are well described by only a single parameter, Sec. 10 discusses a 
recent conjecture stating that the single-parameter liquids are precisely those that 
obey density scaling for the results of high-pressure experiments. Finally, Sec. 11 
gives a brief summary. 



2. The conventional wisdom: One parameter is seldom enough 

The standard "order" parameter theory of glass science was developed by Davies and 
Jones in the 1950's [7l [8]. This theory idealizes the glass transition and treats it 
as a genuine phase transition. In the liquid the "order" parameters are functions 
of pressure and temperature, whereas they are frozen in the glass phase. If Ac p 
is the difference between liquid and glass isobaric specific heat per unit volume at 
the glass transition temperature T g , Akt the liquid-glass difference of isothermal 
compressibilities, and Aa p the liquid-glass difference of isobaric thermal expansion 
coefficients, the Prigogine-Defay ratio II is defined [HI [9] by 

n = AcpAKT 2 . (i) 

Within the Davies- Jones framework one can prove [8) that II > 1 , an inequality 
that has been confirmed in many experiments on quite diverse glass-forming liquids 
[ini El HH H3] • If there is just a single "order" parameter, one has II = 1. Although 
there are glass-forming polymers where II = 1 within experimental uncertainty [14U15) . 
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the vast majority, if not all, glass- forming liquids have II > 1 (typically: 2 < II < 5) [2J. 
If simple first-order dynamics are adopted, the case of a single parameter implies an 
exponential decay towards equilibrium after external disturbances [U [21 [SI HI El [HI El [H] - 
This is rarely observed, a fact that traditionally was seen as a confirmation of the 
conventional wisdom that more than one parameter is required. 

A further classical argument for one parameter not being enough is the well- 
known fact that glass properties are not uniquely defined by, e.g., the density, as one 
would expect if there is just one parameter [3] Qj]. For instance, one can prepare 
glasses with same index of refraction, but different electrical conductivity. This point 
was beautifully illustrated in Kovacs' classical cross-over experiments pTlfTp] . 

In summary: The observed Prigogine-Defay ratios are almost always significantly 
larger than unity, relaxations are almost always noncxponential, and glass properties 
are not just a function of density. This altogether makes a convincing case for there 
generally being a need for more than one parameter, a conclusion that also appears 
natural given the complexities of glass-forming liquids and glass structure. Based on 
this, with few exceptions (e.g., Ref. [TTJ [18l [19l [20]), the matter has not been actively 
discussed for long time. 

3. Questioning the conventional wisdom 

The first point to be noted is that the question of one or more "order" parameters 
is not really well defined in the classical approach, because the glass transition is 
not a phase transition. The fact that the glass transition is a dynamic phenomenon 
- a gradual falling-out-of-equilibrium that inevitably takes place whenever inherent 
relaxation times become longer than experimental times - is well known and well 
understood. This weakens the classical theory where one regards the glass transition 
as a freezing- in process taking place at a particular temperature [21] . 

A related conceptual problem is that the II of Eq. ([1]) is not strictly well defined. 
The changes in specific heat, etc, from liquid to glass are not well defined because of 
two facts: 1) These changes are found by extrapolating the liquid and glass properties, 
respectively, to the transition region. The glass transition temperature, however, is 
not strictly well defined because the glass transition is not a phase transition. 2) The 
glass phase is not well defined - and it relaxes continuously - in principle making 
any measured property in the glass phase a function of time. Many researchers would 
argue that, while this is correct in principle, these effects are minor and not sufficiently 
important to reduce the observed Prigogine-Defay ratios to unity. We take a more 
purist viewpoint, however, and believe that concepts that are not well defined should 
be avoided in a scientific description. 

Recent experiments monitoring ageing of a glass at temperatures around T g 
indicate that in some cases the deviations from equilibrium may be quantified in 
terms of a single parameter. One example is a study where the characteristics of the 
dielectric Johari-Goldstein beta loss peak were used to monitor structural relaxations 
taking place on the alpha-time scale [22J, [23] . To keep things simple only beta loss 
peak frequency and beta maximum loss were monitored, thus providing two numbers 
that depend on structure and temperature. For both sorbitol [22] and tripropylene 
glycol [23j it was found that at any given temperature these two numbers correlate 
linearly. Thus even after a complex thermal history, when returning back to some 
given temperature, beta loss-peak frequency and loss maximum always lie on a line 
characterizing that temperature. An example of this is provided in Fig. 1 showing 
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Figure 1. Temperature-jump experiment for tripropylene glycol monitored via 
the beta loss-peak frequency and loss peak maximum |23| . These two numbers 
depend on structure, and their relaxation monitors structural relaxation (known 
to take place on the alpha-time scale). Starting at 185.0 K temperature was 
first lowered to 183.0 K and kept there for 84 h. Then temperature was changed 
to 181.0 K where it was kept for 140 h. Thereafter temperature was changed 
back to 183.0 K and kept there for another 140 h. The fact that the beta loss- 
peak frequency and loss maximum correlate at any given temperature, also after 
jumping to 181.0 K and back go 183.0 K, indicates that the structure controling 
the beta relaxation may be described by a single parameter. 



loss-peak frequency and maximum loss for the Johari-Goldstein beta process of tri- 
propylene glycol during a temperature cycling around T g . If the structure were 
characterized by more than one order parameter, there is no reason why such a 
correlation should hold. On the other hand, if structure is characterized by a single 
parameter, at any given temperature the two quantities must correlate, and for fairly 
small deviations from equilibrium this correlation would appear approximately linear. 

Other dielectric experiments also indicate that a single structural parameter may 
be sufficient in some cases. Thus studying the shape of the alpha loss peak as quantified 
by the exponent of the best-fit stretched exponential function, it has been shown for a 
number of liquids [24] that when both temperature and pressure are varied, the shape 
depends only of the loss-peak frequency. A simple explanation of this would be that 
there is just a single structural parameter, because if that were the case, this parameter 
would determine both loss-peak frequency and loss-peak shape and consequently these 
two quantities would automatically correlate. 

Richert and Weinstein from a study of the nonlinear dielectric response on glycerol 
showed that, although the dielectric and thermal relaxation times vary throughout the 
liquid, they are locally closely correlated [25] ■ Again, if there is just one parameter 
determining all properties, one would expect that this parameter may fluctuate in 
space, but locally determine both dielectric and thermal relaxation time. 

More evidence comes from computer simulations. In a model of ortho-terphenyl 
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Mossa and Sciortino [26| studied ageing for fairly small temperature steps that were, 
however, large enough to be well outside the linear regime. The simulations showed 
that in configuration space the location of the aging system can be traced back to 
equilibrium states. The authors summarized their findings by stating that for non- 
linear relaxations close to equilibrium "a thermodynamic description based on one 
additional parameter can be provided." 

There seems to be a general understanding in the glass community that there is 
one "order" parameter if and only if the Prigogine-Defay ratio is unity, which happens 
if and only if relaxations are simple exponentials. This is not correct, however, and 
not what one finds from reading the classical papers carefully. Goldstein in his 1964 
review, for instance, noted that in some situations with several parameters that are 
mathematically constrained, the Prigogine-Defay ratio may be unity. In this situation 
"it is really a matter of taste" [27] whether one prefers to speak of many (constrained) 
parameters or of a single generally non- exponential "order" parameter. Thus the 
observation of non-exponential relaxations does not imply that there must be more 
than one parameter. 

In 2006 Schmelzer and Gutzow revisited the Prigogine-Defay ratio and the 
question of the number of "order" parameters [20]. Assuming just a single parameter 
the dynamics of which follow the classical framework of the thermodynamics of 
irreversible processes, they showed that the standard Prigogine-Defay ratio obtained 
by extrapolating from glass and liquid to T g must be larger than unity. This result 
seriously questions the prevailing understanding of the glass community referred above, 
and it emphasizes the need for further work. 

The first one should do in reconsidering the "order" parameter question is to 
make sure that the problem is well defined. As mentioned, the standard Prigogine- 
Defay ratio II of Eq. ([]]) is not well defined. As became clear in the 1970's 
|H HI, it is possible to define a version of II that is well defined. This is 
done by referring exclusively to properties of the equilibrium viscous liquid phase 
and its linear responses. In this phase thermodynamic properties are generally 
frequency dependent, and the high-frequency limits correspond to glassy behavior 
where structural relaxations do not take place. If c p (w) is the frequency-dependent 
isobaric specific heat per unit volume [28], etc, this leads to the following rigorous 
definition of the Prigogine-Defay ratio for the metastable equilibrium viscous liquid 
at any temperature T: 

n _ [c p (oJ -> 0) - c p (uj -» oo)] [k t (u> -* 0) - k t (uj -» oo)] ^ 
T\a p (uj — > 0) — a p (uj — > oo)] 2 

4. Linear thermoviscoelasticity 

From now on we turn the focus exclusively to the metastable liquid phase with 
no reference to the glass phase. This limits the discussion compared to what is 
standard in glass science, but has the advantage of making all concepts rigorously 
well defined. Linear thermoviscoelasticity deals with the frequency-dependence of 
thermodynamic properties and their coupling to frequency-dependent mechanical 
properties. It is understood that, in principle, only infinitesimal perturbations are 
applied, thus ensuring linearity. In the simplest (isotropic) theory there are two 
fundamental "energy bonds," a thermal and a mechanical. An energy bond has an 
"effort" variable and a "displacement" variable 29, 30, 31j. The thermal energy bond 




Figure 2. The two fundamental energy bonds for a system described by standard 
thermodynamics. One energy bond is thermal; here the "effort" is temperature T 
and the displacement variable is the entropy S; thus if S = dS/dt is the entropy 
flux into the system, the rate of energy transferred into the system is TS. The 
second energy bond is mechanical; here the effort is negative pressure —p and the 
displacement variable is the volume V; if V is the volume flux into the system, 
the rate of energy transferred into the system is —pV. 



is characterized by entropy S as the displacement variable and temperature T as 
the effort, for the mechanical energy bond the displacement variable is the volume 
V and the effort is the negative pressure, —p. The product of the effort and the 
differential displacement variable gives the energy transferred into the system from 
its surroundings. Thus the two energy bonds (Fig. 2) simply express the well-known 
fundamental identity dE = TdS — pdV . 

For infinitesimal perturbations around equilibrium with angular frequency to, if 
one imagines controlling the effort variables and measuring displacement changes, 
and if the usual complex notation is adopted where, e.g., T(t) = To + 5T(t) with 
ST(t) = He[6T exp(iu)t)], linearity is expressed in the following relation where ds is 
entropy change per unit volume and Sv is relative volume change: 



/ Ss(u;) \ _ ( c p (u;)/T a» \ / ST(u) \ , . 

{ 5v(oj) ) ~ \ a p (uj) k t (u) J { -5p(u>) J ' (6) 

The response matrix is sometimes termed the thermal compliance matrix. Its 
symmetry expresses Onsager reciprocity, reflecting the fundamental fact that there is 
time reversibility on the microscopic level [32 [33] . 



5. The completely general case: Three independent thermoviscoelastic 
response functions 

How many independent thermoviscoelastic response functions exist? From the four 
variables, entropy, temperature, volume and pressure, one may choose any two as 
"control" variables. Usually, one chooses one control variable from each energy bond. 



Glass-forming liquids: One or more "order" parameters? 



7 



There are thus four natural choices of control ("input") variables, the two remaining 
are the measured ( "output" ) variables. For each choice there is one response matrix as 
in Eq. © 34 . These four matrices are all symmetric by Onsager reciprocity, leaving 
12 frequency-dependent response functions. These are not independent, however; if 
one matrix is known, the three others are easily calculated from it by isolating the 
output variables in question on the left hand sides of two equations. Thus there are 
only the three independent response functions, for instance those of Eq. ([3]). 

The above statement is true in complete generality. Papers by Moynihan and 
others of the 1970's showed, however, that in the "order" parameter description there 
are really only two independent response functions [H 01 [35] (see also Ref. [36]). 
The formula for calculating the third response function in terms of the two others 
involves analytic continuation (next section). The situation is analogous to that of 
the Kramers-Kronig relation which allows one to calculate the imaginary part of a 
response function in terms of its real part, but only if the latter is known at all 
frequencies. 

6. The general stochastic case: Two independent thermoviscoelastic 
response functions 

In this section we summarize the master equation description of viscous liquid 
dynamics [34] [3TJ [38] and show that it implies that there are just two independent 
response functions. More precisely, it is shown that knowledge of c p (u>) and a p (u>) at 
all frequencies allows one to calculate Kt{u) except for an overall additive constant 
giving the high-frequency limit. This is equivalent to the above-mentioned result 
derived long ago [TJ [3] [35], but now in a setting that is explicitly consistent with 
statistical mechanics. 

In a master equation there are states and stochastic transitions between the states. 
A complete description is provided by the set of probabilities {Pn} that the system is in 
state n. This is an ensemble description making it possible to calculate all properties, 
including the entropy. Following Ref. [34] we shall think of each state as an inherent 
state in the sense of Stillinger and Weber [39] (i.e., a potential energy minimum in 
configuration space), but other state interpretations are also possible. Each state has 
the vibrational Gibbs free energy G n (T,p). The ensemble Gibbs free energy that 
includes the probability dependence is given [M] [35] by 

G{T,p,{P n }) = Y d P n (G n {T,p) + fc s TlnP„) . (4) 

n 

From this one finds the ensemble volume and entropy by the usual thermodynamic 
relations V = dG/dp and S = -dG/dT. 

The master equation dynamics are given by first order equations in time that are 
mathematically similar to those of the classical "order" parameter description of glass 
science: 

Pn = ^W nm P m . ( 5 ) 

m 

The main difference to the "order" parameter description is the constraint P n = 1 
and that the present formalism ensures consistency with statistical mechanics. 

The rate matrix W depends on T and p and changes slightly when these 
variables are perturbed by small time-dependent variations. The same applies for 
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the equilibrium probabilities, P£ q oc exp[— G n (T, p)/ksT]. According to the principle 
of detailed balance, which ensures consistency with statistical mechanics as well as 
time-reversal invariance, the equilibrium probabilities [321 133j obey 

W nm (T,p)P%(T,p) = W mn (T,p)P^(T,p). (6) 

Here temperature and pressure may be arbitrary functions of time. For periodic 
infinitesimal perturbations from equilibrium the dynamics are perturbed via the 
transition matrix's dependence on pressure and temperature. The equilibrium 
probabilities at p — Pq and T — T are denoted by P° and the transition matrix 
at this state point is denoted by W°. 

If Qi is entropy and Qi volume, solving the resulting system of equations leads [34 
to the following expression for the compliance matrix of Eq. §3§ (where and ^p^- 
are evaluated at (Tb,po), the matrix A(ui) is defined by A{lo) = (W° — iu>)~ 1 W°P° , 
a,0=l,2): 



Note that A(lo) — > for uj — > oo; thus J a p{Lu) — > J^g for uj — > oo. Introducing 

the matrix Y nm = (P°) _ ^W^ TO (P^ l )i , the detailed balance requirement Eq. ([6]) 
implies that Y is symmetric. In terms of Y, the matrix A(u>) is given [34] by 
A(to) — RYiY — iu)^ 1 R where R nm = {Pn)^S nm . Thus for the relaxing part of 
the compliance matrix A J = J — J°°, if dQ a is the vector whose n'th component is 
dQ a /dP n (T Q , PQ ), one has [34] 



AJ af} (u) = (RdQ 



Y 



Y-ico 



RdQ a) . (8) 



Adopting the standard "ergodicity" assumption that all states are connected by 
some path of intermediate states, the matrix Y has a one-dimensional eigenspace 
corresponding to the eigenvalue zero whereas all other eigenvalues are negative [3^1l33j . 
If the eigenvectors of Y corresponding to all the negative eigenvalues are denoted by 
\ipj) with corresponding eigenvalue — 1/tj, Eq. (jSj) implies 

aj q/3 h = ^(i?ag a \^)(^\RdQf,)—^— . (9) 

Since (— 1/t,-)/(— 1/tj —iu) — 1/ '(1 + iwTj), changing to a continuous notation Eq. © 
becomes 

AJ,*(a,)= r^aMrfr, (10) 

7 1 + IWT 

where the functions g a (r) are real, but not necessarily positive. 

By reference to the theory of analytic functions we show below that not all 
three functions of the compliance matrix are independent. This is intuitively obvious 
already from the fact that the three compliance functions are determined by the two 
functions gi(r) and g% (r). More precisely the argument goes as follows. The three 
compliance functions AJ a ^((j) are analytic in the complex positive half plane, i.e., 
where Re(w) > 0. Knowledge of such a function at all real, positive frequencies by 
analytic continuation uniquely determines the function in the complex plane. Equation 
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(|10p shows that there is a branch cut along the positive imaginary frequency axis. 
Given that AJ a p(u) — > for uj — > oo, the pole distribution on the branch cut uniquely 
determines the compliance function. More specifically, Eq. (fTI)]) implies that 

AJ 22 H = ^ lim {(i + iw V)^M\ 1_ dr. (11) 

J uj>->i/t I £\J\i(uj') J 1 + IUJT 

Thus knowledge of Ac p (w) and Aa p (w) at all real, positive frequencies implies 
knowledge of Akt(lo). Similarly, knowledge of Akt(^) and Aa p (w) at all real, positive 
frequencies implies knowledge of Ac p (u>). 

7. The single-parameter case: One independent thermoviscoelastic 
response function 

The compliance matrix J a p(uj) reflects both the relaxing responses (completely 
characterized by AJ(w)) and the instantaneous responses given by the high-frequency 
limits. Switching to the time domain, if the relaxing responses of the two energy bonds 
are always proportional, i.e., controlled by a common variable Se(t), the entropy and 
volume responses per unit volume are given by expressions of the form 

Ss(t) = 7l fe(t) + J£6T(t) - J£6p(t) 

Sv(t) = 72 Se(t) + J%6T(t) - J%6p(t) . (12) 

We refer to this situation as that of a single "order" parameter [33] and proceed 
to show following Ref. [33] that in this case there is basically just one compliance 
function. Note that no reference is made to the properties of the glassy state. 
For periodically varying fields Eq. (p~2|) implies 

8s(u) = 7i<5eH + J™8T{u)) - J^5p{u) 

Sv(u) = 72<5eH + J£6T(lu) - J^5p{u) . (13) 

The e-parameter may be expanded to first order as follows: 

5e(lu) = Ai(w)5T(w) - A 2 (uj)6 P (lu) . (14) 

Substituting Eq. (|14[) into Eq. (|13|) and using the symmetry of the compliance matrix 
leads to the identity 7iA 2 (w) + J£? = 7 2Ax(w) + J£f. For the imaginary parts this 
implies 

A]>) = A|H^ (15) 
7i 72 

When two analytical functions both with branch cuts on the positive imaginary axis of 
the complex w-plane have same imaginary part, they are identical except for an overall 
additive constant. The latter is zero, because the fact that the two functions give the 
relaxing part of the responses implies that they both go to zero for w — > oo. Thus 
h-xipS) oc A2(u>). By considering the constant pressure and constant temperature cases 
it now follows easily from Eqs. (flU]) and (Tj"4"|) that AJn(w) oc AJi2(w) oc AJ22(w), or: 

Ac p {uj) oc Actp(u>) oc Akp(w) . (16) 

In conclusion, in the case of a single "order" parameter (Eq. (TT2"|) there is basically 
just one independent thermoviscoelastic response function, i.e., knowledge of one of 
them implies knowledge of the two others except for the overall additive constant 
giving their high-frequency limits. 
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8. "Dynamic" Prigogine-Defay ratio: A single-parameter test 

In principle, in order to test experimentally whether or not a single "order" parameter 
suffices, one measures the three response functions of the compliance matrix to test 
whether the relaxing parts are proportional (Eq. P6[) ). This, however, requires wide- 
frequency measurements of the thermoviscoelastic response functions, and there are 
yet no measurements of all three thermoviscoelastic response functions on a glass- 
forming liquid. (Even the isobaric frequency-dependent specific heat c p (u>) has not 
yet been measured reliably [28j . The problem is that, because frozen- in stresses relax 
on the same time scale that the enthalpy relaxes; establishing truly isobaric conditions 
is difficult and in most experimental setups the stress tensor is not diagonal.) 

Even when methods have been developed for measuring the compliance matrix 
of Eq. ([3]) , one may still expect that initial measurements cover only a rather limited 
dynamic range. This leads to the question: Is it still possible to test the single- 
parameter conjecture Eq. (fT"2f ? This question was discussed in a recent publication 
[34] where it was shown that, in fact, measurements at one single frequency are enough 
to test the single-parameter conjecture. Of course, one can never prove that a single- 
parameter description is correct in an absolute sense - it is all a matter of investigating 
how good such a description is. In the above-mentioned recent paper [34] it was shown 
that a "dynamic" Prigogine-Defay ratio At p (w) > 1 exists with the property that, if 
this quantity is unity at one frequency, it is unity at all frequencies - which happens if 
and only if a single-parameter description applies. The dynamic Prigogine-Defay ratio 
is given by the imaginary parts of the three thermoviscoelastic response functions [34] 
as follows: 

In order to minimize uncertainties measurements should preferably be taken at a 
frequency around the alpha loss-peak frequency, because only here the imaginary 
parts are significantly different from zero. We expect that if At p {w) is close to unity 
in the main relaxation region (e.g., below 1.1), a single-parameter description applies 
to a good approximation. 



9. Results from computer simulations 



Recently thermal equilibrium fluctuations were studied in computer simulations of 
various liquids [401 SI] ■ I n many cases it was found that in constant temperature and 
volume simulations (the so-called NVT ensemble) pressure and energy fluctuations 
correlate strongly. More accurately, this applies for the configurational parts of 
pressure and energy, the "virial" and the potential energy. (The kinetic parts of 
pressure and energy - the ideal gas pressure at the given density and temperature, 
and the kinetic energy - trivially correlate 100%, but with a different proportionality 
constant.) As an example, Fig. 3 shows the thermal fluctuations of virial and potential 
energy for a standard Lennard-Jones liquid. 

Liquids for which these quantities correlate strongly in their fluctuations are well 
described by a single order parameter [33]. Intuitively this may be understood by 
reference to Eq. (fT2|) considered without perturbations (6T(t) — Sp(t) = 0) which, if 
assumed to describe also the fluctuations, shows that entropy and volume fluctuations 
are 100% correlated. Thus one expects that the dynamic Prigogine-Defay ratio is close 
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Figure 3. Thermal equilibrium fluctuations of potential energy and viral (the 
configurational part of pressure) for a standard Lennard- Jones liquid [41] , The 
fact that these two quantities correlate strongly shows that, as regards the 
configurational degrees of freedom, a single-parameter description is quite good 
for the thermoviscoelastic behavior. For highly viscous liquids the time scale 
separation between the slow configurational degrees of freedom and the remaining 
implies that these correlations (that we have also seen, e.g., in simulations of the 
highly viscous Kob-Andersen binary Lennard-Jones mixture) implies that the 
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Figure 4. Correlation coefficient between virial (volume times the configurational 
part of the pressure) and potential energy thermal equilibrium fluctuations for a 
number of liquids evaluated by computer simulations 141 1 . The liquids represented 
are: LJ: Standard Lennard-Jones, Exp: Monatomic liquid with exponential 
repulsive forces, Dumb-bell: A molecule model of two atoms of unlike size, BLJ: 
The Kob-Andersen binary Lennard-Jones liquid, Methanol, and SPC/E Water. 
The last two are hydrogen bonding and do not show significant correlations; the 
other liquids do. It has been argued that virial and potential energy give the 
slowly fluctuating parts of the pressure and energy |41l ; thus whenever the former 
quantities correlate strongly, to a good approximation the liquid may be regarded 
as described by a single parameter. 
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Figure 5. Correlation coefficients for a number of glass-forming liquids between 
virial and potential energy thermal equilibrium fluctuations as function of 
temperature (in reduced units). The liquids represented are: LJ: Standard 
Lcnnard- Jones, Exp: Monatomic liquid with exponential repulsive forces, Dumb- 
bell: A molecule model of two atoms of unlike size, BLJ: The Kob-Andersen 
binary Lennard- Jones liquid, Methanol, and SPC/E Water. The figure shows 
the same systems as those of Fig. 4 studied by computer simulations. The 
two hydrogen-bonding liquids, water and methanol, show poor correlation, the 
remaining systems are all strongly correlating. For the latter the correlation even 
increases as temperature rises; this is because at high temperature the particles 
approach each other more in collissions than at low temperatures and the inverse 
power law description of the works better the closer the particles are. 



to unity for such "strongly correlating liquids." Figure 4 shows that this is the case for 
the Lennard-Jones liquid as well as for a number of other glass-forming liquids. Water 
and methanol are interesting exceptions that do not show strong correlations between 
virial and potential energy fluctuations (Fig. 5); thus for these two hydrogen-bonding 
liquids a single-parameter description does not apply. 

10. A conjecture: Strongly correlating liquids obey density scaling and vice versa 

The last five years large amounts of data on the behavior of glass- forming liquids under 
pressure have been published. The motivation is that by not just varying temperature, 
but pressure as well, much more information may be learned about these systems. 
Generally, the liquid relaxation time r, which is basically the Maxwell relaxation time 
or the alpha loss peak frequency, depends strongly on both temperature and pressure, 
increasing with lowering temperature or raised pressure. This is not surprising. A new 
and significant finding 031 03] , however, is that if p is the density, many liquids 
obey "thermodynamic" or "density" scaling, i.e., the function r(T,p) may be written 




Both the function F and the exponent x depend on the liquid in question. This 
expression has mainly been tested on glass-forming molecular liquids, the systems 
that are most easily accessible. For hydrogen-bonding liquids like glycerol or sorbitol 
the x's initially reported were anomalously small [44], but it now appears that the 
reason is that density scaling does not work very well for hydrogen-bonding liquids 
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Recently, Coslovich and Roland presented computer simulations of binary 
Lennard- Jones type systems where, however, the exponent of the repulsive term of 
the potential was varied, taking the values 8, 12, 24, and 36 [35]. Such system may 
be cooled to low temperatures where the viscosity is very large, without crystallizing. 
Their simulation results obey the density scaling expression Eq. (fl8)) . which by itself 
is an interesting finding. Even more interesting is the fact that the exponent x may be 
related to the effective exponent describing the approximate power law of the potential. 
For the standard binary Lennard-Jones case, for instance, this exponent is not 12 as 
naively expected, but a number close to 18 depending on the precise choice of fitting 
criteria [4T] , 

In Ref. [41] some of the present authors previously found that there are strong 
energy-pressure correlations whenever the repulsive part of the interaction is well 
described by an inverse power law. Since this seems also to be the criterion for a liquid 
obeying density scaling (Eq. (fig)) ), an obvious conjecture is [H] that: A glass-forming 
liquid is strongly correlating if and only if it obeys density scaling. Two liquids that 
in computer simulations were not strongly correlating are water and methanol [41 J , 
and we surmise that hydrogen-bonding liquids generally are not strongly correlating. 
The argument is that the existence of "competing interactions" (van der Waals forces 
as well as the directional hydrogen bonds) destroy significant correlations, implying 
that hydrogen-bonding liquids are not well described by a single "order" parameter. 
This is consistent with the finding that hydrogen-bonding liquids do not obey density 
scaling. 

If this conjecture is correct, by virtue of their simplicity the class of strongly 
correlating liquids provides an obvious starting point for theories for viscous liquids 
and glass formation. It would be obvious to further conjecture that also covalently 
bonding liquids are not strongly correlating, again due to the directional nature of 
the bonds. Many of these systems have fairly low fragility. Low-fragility liquids 
are traditionally thought to be simple (e.g., have to almost exponential relaxations 
if the liquid is almost Arrhenius). We here conjecture almost the opposite, namely 
that many high-fragility liquids in a certain sense are simpler than many low-fragility 
liquids. 

11. Summary and final remarks 

We have argued that the old discussion of one or more "order" parameters deserves 
to be revitalized. There are indications that at least some glass formers may be well 
described by a single order parameter as regards their linear thermoviscoelasticity. 
It is important to emphasize that no claim is made that the molecular structure is 
completely characterized by a single number. We now have an experimentally useful 
criterion for whether or not a single-parameter description is accurate. Computer 
simulations confirm that some model liquids are well described by a single parameter; 
these liquids are referred to as "strongly correlating." Since hydrogen-bonding liquids 
do not show these correlations, we expect that liquids with directional bonding are 
not well described by a single parameter, whereas van der Waals bonded liquids are. 
Thus we conjecture that for van der Waals liquids the relaxing parts of the three 
thermoviscoelastic response functions of Eq. ([3]) are all proportional, whereas for 
hydrogen-bonding liquids this is conjectured not to be the case. This prediction can 
be tested once methods have been developed to measure the full thermoviscoelastic 
compliance matrix. 
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